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Abstract. A binary compact object early in its inspiral phase will be picked up by 
its nearly monochromatic gravitational radiation by LISA. But even this innocuous 
appearing candidate poses interesting detection challenges. The data that will be 
scanned for such sources will be a set of three functions of LISA's twelve data streams 
obtained through time-delay interferometry, which is necessary to cancel the noise 
contributions from laser-frequency fluctuations and optical-bench motions to these data 
streams. We call these three functions pseudo-detectors. The sensitivity of any pseudo- 
detector to a given sky position is a function of LISA's orbital position. Moreover, at a 
given point in LISA's orbit, each pseudo-detector has a different sensitivity to the same 
sky position. In this work, we obtain the optimal statistic for detecting gravitational 
wave signals, such as from compact binaries early in their inspiral stage, in LISA data. 
We also present how the sensitivity of LISA, defined by this optimal statistic, varies 
as a function of sky position and LISA's orbital location. Finally, we show how a real- 
time search for inspiral signals can be implemented on the LISA data by constructing 
a bank of templates in the sky positions. 
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1. Introduction 

The commissioning of Earth-based long-baseline gravitational wave (GW) interferom- 
eters has finally come to fruition a little less than three decades since the discovery of 
the Hulse and Taylor binary pulsar in 1974 pQ and the subsequent confirmation of the 
emission of gravitational waves (GWs) by that system. [2] On the other hand, the Laser 
Interferometric Space Antenna (LISA) is being designed as a space-based detector to 
observe low-frequency GWs (in the milli- to deci-hertz band) to complement the high- 
frequency observations (in the deca- to kilo-hertz band) of its Earth-based counterparts. 
Many of the sources of low-frequency GWs fall in the category of compact binary ob- 
jects, which include white dwarfs, neutron stars, and black holes. During the early stages 
of the inspiral, these compact binaries will emit almost monochromatic low-frequency 
gravitational waves, which will produce detectable signals in the LISA data. 

LISA will comprise three spacecrafts located at the vertices of a nearly equilateral 
triangle (as shown in Fig. with the side lengths equal to 5 million kilometers. [3 3 Each 
craft will house a couple of laser-mounted optical benches and proof masses and will 
freely fall around the Sun in an orbit that lies on a plane slightly tilted with respect to 
the ecliptic. The tilt will vary from one craft to the other. By inclining the plane of the 
triangle to the ecliptic by a constant angle of 60°, the side- length of the triangle can be 
maintained at a fixed value. In the process, the triangle will complete one spin about 
its normal while its centroid, trailing behind the Earth by about 8.3 million kilometers, 
completes one orbit around the Sun. 

The effect of a GW on LISA will be to change the physical distances between 
its freely falling proof masses. This change will be registered as fractional frequency 
shifts in the six laser beams exchanged among the three space-crafts. The shifts 
will have additional contributions from various noise sources. These include two 
primary ones, namely, the laser-frequency fluctuations (contributing a fractional shift 
of about 10 -13 Hz _1//2 ) and the optical-bench motions (with a fractional shift of about 
1(T 16 Hz~ 1/2 ). In order to detect the GWs (with a strain around 1(T 21 Hz~ 1/2 ), it 
is imperative that these noises be mitigated by several orders of magnitude. A data 
analysis technique for achieving this goal was accomplished by Armstrong, Estabrook, 
and Tinto.|l]. They showed that by combining appropriately the time-delayed versions 
of these six data streams, with six additional ones arising from laser beam exchanges 
between adjacent optical benches on each craft, one can eliminate the two primary 
noises, thus, rendering the LISA data analyzable for GW signals. 

In this paper, we formulate a strategy for detecting nearly monochromatic 
gravitational waves from inspiraling compact binary objects in the LISA data. There 
are two complementary aspects to such a strategy. The first is to deduce the maximum 
number of noise-independent detectors that LISA offers. And the second is to construct 
the appearance of a GW signal in them. This allows one to match the data from 
these detectors with a template of the expected signal in them. Whether a match is 
strong enough to warrant a detection is then decided based on the rate of false alarms 
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at the level of that match. For the problem of detecting low-mass compact binaries, 
involving white-dwarfs and neutron stars, the waveforms obtained by the quadrupole 
approximation [Sj suffice for accurately modeling the expected signal. This waveform 
allows for a slightly "chirping" source, i.e., a source whose orbital frequency Qq is 
increasing at a rate Qq -C Qq/T, where T is the observation period. We will use this 
waveform clS db template in our detection strategy. 

The layout of the paper is as follows. In Sec. El we enunciate the three noise- 
independent data combinations, or "pseudo-detectors" , that were first obtained in Ref. 
[S| by combining the frequency shifts of the twelve data streams exchanged among the 
LISA spacecrafts through time-delay interferometry (TDI) In Sec. El we obtain 
the form of the gravitational-wave strain caused by a compact binary source in these 
three pseudo-detectors. This provides the templates required to search for signals from 
such sources in the LISA data. We then study the antenna patterns of these pseudo- 
detectors at different points in LISA's orbit, emphasizing how their relative sensitivities 
to a mildly chirping source vary as a function of sky position and source frequency. 
The study bears out the fact that noise-independence of detectors is a property distinct 
from their geometric independence. Indeed, at GW wavelengths larger than about 0.1 
AU, which we will term as the long- wavelength limit, the strain of the third pseudo- 
detector tends to the difference of the strain of the first two. We then derive the optimal 
statistic for detecting gravitational waves from (non-spinning) compact binary inspirals 
by coherently combining the data of the three pseudo-detectors in Sec. EJ 

The statistic obtained here tracks the Doppler modulation of the source frequency 
induced by the motion and time-varying orientation of LISA with respect to that 
source. [0] For Gaussian noise, we derive the probability distribution of our statistic, 
which can be used to compute the signal-to-noise ratio (SNR) of a candidate event. 
We use this distribution to predict the behavior of false-alarm and detection rates as a 
function of the detection threshold set for the statistic. In Sec. El we construct the metric 
on the space of parameters that allows one to estimate the fractional loss of SNR for a 
given mismatch between the template and the signal parameters. Using this metric we 
esimate the number of templates that will be required to search the full astrophysically 
relevant volume of the parameter space while suffering a loss in SNR of no more than 
3%. We also estimate the computational speed required to implement such a search in 
real time. We briefly summarize the results obtained in Sec. El especially, stressing the 
applicability of our formalism to searches in the recently found second-generation TDI 
data combinations. IT2"| HB*] Note that in the expressions appearing in this paper, we 
set the gravitational coupling constant (G) and the speed of light (c) to unity. 

2. The pseudo-detectors 

As illustrated in Fig. [TJ LISA consists of three spacecrafts, labeled i =1, 2, and 3, located 
clockwise at the correspondingly labeled vertices of an almost equilateral triangle. Let 
the arm- lengths of this triangle be L%, L2, and L3, such that Li is the length of the arm 
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Figure 1. LISA consists of 3 spacecrafts located at the vertices of an equilateral 
triangle. These craft exchange six elementary data streams, labeled Ui and — The 
Ui beams propagate clockwise, whereas the — V% propagate counterclockwise. 



facing the vertex i. For specifying the orientation of each arm, we assign unit vectors 
hi along them, such that their directions are oriented anticlockwise about the triangle. 
Each spacecraft will have two optical benches (denoted by % and i* in Fig. |TJ) equipped 
with independent lasers and photodetectors. Thus, each spacecraft will shoot two beams 
towards the other two spacecrafts, respectively, resulting in six one-way beams along 
the 3 arms. 

The effect of an impinging GW is to cause a shift in the frequency of the laser 
beams. But a GW is not the only source causing such a shift. The laser-frequency 
fluctuation is another source. If vq is the central frequency of all the lasers in LISA, 
then the fractional shift caused by such fluctuations in the beam originating at optical 
bench i is 

cm s *a®. . (1) 

Similarly for the beam from the bench i*. What is measured, however, is the frequency 
fluctuation in the beam from one bench relative to that in the beam from a bench in 
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one of the two other vertices. This way one obtains three "clockwise" oriented streams, 
Ui, and three "anticlockwise" oriented streams, —V^ The frequency fluctuation in the 
beam from bench {% — 1) relative to that in the beam originating at bench i* is termed 
Ui. \ Therefore, U\{t) = C 3 (t — L 2 ) — Ci(t) is the data stream measured by beating the 
beam transmitted by bench 3 against that of bench 1*, measured at time t at bench 1. 
The remaining two streams, namely, Ui and U 3 , can be obtained by cyclic permutation 
of the indices in the U\ expression. (Thus, U 2 {t) = Ci(t — L 3 ) — C 2 {t), and so on.) Three 
more data streams, termed — Vi, are obtained by measuring the frequency fluctuation 
in the beam from bench (i + 1)* relative to that originating at bench i. For instance, 
—Vi(t) = C 2 (t — L 3 ) — Ci(t). Once again, the remaining two — Vi streams can be 
obtained by cyclic permutation of the indices in the above expression for —V\. (Hence, 
-V 2 {t) = C 3 (t - Li) - C 2 (t), and so on.) 

The fractional shift C 3 (t — L 2 ) is constructed from C 3 (t) by shifting back the latter 
stream in time by an amount L 2 . For brevity of expressions, we introduce the time-shift 
operator Q by its action on a data stream x{t) as in: 

&x(t) = x(t - Li) , (2) 

where the label i denotes the arm along which the time-shift is affected. One can thus 
define the 6 inter-craft streams as follows |T4*]: 

U\ = C2C3 — Ci j U 2 = C3C1 — C 2 , U 3 = C1C2 — C3 j 
V 1 =C 1 -£ 3 C 2 , V 2 = C 2 -( 1 C 3 , V 3 = C 3 -( 2 C 1 . (3) 

Note that the effect of such a shift on the Fourier components, x(f) = xii^e^^dt, 
of the data is to change them to e~ 2m ^ Li x(f), where / is a frequency variable. Therefore, 
the effect of the operator Ylt=i C™' 011 a data stream is to change its Fourier transform 
by the factor e~ 2m ^^ iaiLi . 

These data streams, however, are expected to suffer from several noise sources, two 
of which, viz., the laser- frequency fluctuations and the optical-bench motions, tower over 
the others. The two other noise sources are the photon-shot noise and the fluctuations 
in the motion of the proof masses. The laser-frequency fluctuations directly influence 
the fractional frequency-shifts, Ci, and, therefore, the data streams U and V.. The 
effect of the optical bench motions on the Ci is additive: Let the random velocities of 
the optical benches be Vj and Vj. . This causes a Doppler shift in the frequency of the 
lasers mounted on those benches, which in turn modifies the Ci to 

Ci Ci - n 3 • vi , 

Ci* Ci* + n 2 • vi* . (4) 

The remaining Ci and Ci* can be inferred by cyclically permuting the indices in the 
above expression. If u ifi * are the random velocities of the proof masses on benches i and 

X Note that the indices i and i± 1 can take only 1, 2, and 3 as values. These three numbers are ordered 
clockwise in Fig. ^ By convention, whereas i + 1 equals the number next to i while going clockwise in 
that figure, i — 1 equals the number preceding i. E.g., when i = 3, we take i — 1 = 2 and i + 1 = 1; 
when i = 1, we take i — 1 = 3 and i + 1 = 2. 
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i*, then it is easy to see that the value of U\ and V\ gets affected by additional terms 
2n 2 ■ u r * and 2n 3 • u 1; respectively. Finally, when there is a GW signal present, these 
streams will receive additional contribution owing to the fractional frequency shifts 
caused by it. Thus, in the presence of the above noise sources and signal, the data 
streams get modified as 

Ux - U x = C 2 (C 3 - n 2 • v 3 ) - (tfi. + n 2 • vi.) + 2n 2 • Ul * + < ot + U? w , 
Vt V x = -C 3 (C 2 . + n 3 ■ va.) + (d - n 3 • Vl ) + 2n 3 ■ u a - < ot - 1^5) 

where n^°y. are the photon-shot noises and Uf W and are the GW signals present 
in the Ui and Vi data streams, respectively. The four remaining [/, and Vi streams can 
be obtained from the above expressions by cyclic permutations of the indices. In the 
rest of the paper the Ui and Vi will refer to these modified streams. When a (GW) 
signal is absent, they will still be given by these modified expression, but with Up w = 
and V i GW = 0. The form of U GW and V^ GW in the presence of a signal will be explored 
in the next section. 

In addition to the six inter-craft data streams discussed above there is 
supplementary information available about the noise sources in the intra-craft beams 
exchanged through the optical fibers connecting two adjacent optical benches % and i* 
located in the zth craft. In all there are six intra-craft beams, two per bench pair. But 
the two intra-craft beams within a craft can be beaten against each other to produce 
a single stream that is directly relevant to noise suppression. To wit, at craft 1, by 
beating the frequency of the beam from bench 1 relative to that on bench 1*, one forms 
the stream, 

W x = (Ci - n 3 • vi) - (Ci» +n 2 • Vi*) + n 3 • u x + n 2 ■ ui» . (6) 

Two other intra-craft data combinations, W2 and W3, can be obtained by the cyclic 
permutation of indices in the above expression. Note that these intra-craft streams will 
bear negligible influence from any impinging gravitational wave. Nevertheless, as we 
explain below, they offer information on laser-frequency fluctuations and optical-bench 
motions that can be used to render the Ui and Vi streams essentially free of any noise 
from these two sources. Together with the U and Vi, the Wi form a a total of nine data 
streams that a data analyst has recourse to in hunting for a GW signal in LISA. 

Following the work of Tinto and Armstrong 7\, it was shown by Dhurandhar et al. 
P^j that by acting on the 6 inter-craft streams, Ui, Vi, and the 3 intra-craft streams, 
Wi, with certain polynomials, pf, qf, and rf, of the time-shift operators, Q, one can 
form several combinations of time-delayed data streams, 

= E \-Pt V M + rfViit) + rfWM] , (7) 

i=l 

that have the laser-frequency and the optical-bench motion noise eliminated. Above, 
A labels the different combinations so obtained. The above technique of constructing 
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such data combinations of pseudo-detectors is called time-delay interferometry. The 
pseudo-detectors can be recast as 

Trace [e A ■ Z] , (8) 



x A 



where 



e 





A =\ q? Q A q A and Z = V 2 U 2 W 2 . (9) 



Thus, for a given choice of the data streams Ui, Vi, Wi (and, therefore, the matrix Z), the 
matrix e A of the time-shift polynomials forms a representation of the pseudo-detectors 



x A . 



Of the several possible pseudo-detectors, only 3 are linearly independent and have 
a non- vanishing GW strain in them, in general [HI ■ The three that will be discussed 
here are defined by their corresponding e A : 

i c 

e 1 = I 1 + 2C 1 - C -2-C I , e 2 = I -1 1 + C -( I ,(10) 





-l + c 2 o 



and 



e 3 = ( 1 1 1 I . (11) 

V -l-C -i-c -i-c 

In the above expressions, it is assumed that all arm lengths are almost identical. 
Therefore, £i ~ £2 — C3 = C It is important to note that these data combinations 
diagonalize their noise covariance matrix |16J and, therefore, are also noise- independent. 

Although the laser- frequency noise is eliminated in the x (t), there is still present 
the noise associated with the acceleration of the proof masses onboard each craft and 
the photon-shot noise. In general, 

x A (t) = n A (t) + h A (t) , (12) 

where 

i=i 

is the gravitational- wave strain in pseudo-detector A and n (t) is the time-delayed sum 
of the remaining noise components. One typically assumes that these components and, 
hence, the total noise, has a Gaussian probability distribution with a zero mean. Their 
variance-covariance matrix elements are given as follows. For the photon-shot noise, 
these elements are 



it 



nT U (f) nf°V) = nt ot *(f) npHf) = n*f *(/) ngf (/') = , (14) 
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for any % and j. Above, p shot (y) is termed as the one-sided power-spectral density 
(PSD) of the photon-shot noise. Furthermore, one assumes the proof-mass noise to be 
isotropic, such that: 

1 



(n, • u*(/)) (n fe • u,(/')) = (fi* • fij.(/)) (A* ■ u,.(/')) = -P proof ,5(/ - f 



(Vu**(/)) (n fc .u i (/0)=(6 4 -uj?(/)) (fi*-u I .(/ , )) = , 

for any z, j, A;, and Z. Also, the covariance of the shot noise in any data stream 
with the noise in the motion of any proof-mass is zero. It is estimated that 
pshot = Lg x io- 37 [//1Hz] 2 Hz _1 and that the proof-mass noise PSD is P proof = 
2.5 x 10~ 48 [//lHz] _2 Hz _1 .l3 j While the proof-mass noise enters the very low frequency 
band of LISA, the shot noise enters the higher end of LISA's sensitivity band. 

It is now possible to deduce the noise PSD, p( A \f), of each of the pseudo-detectors 
from the above expressions. It follows from them that in the absence of a signal, each 
pseudo-detector is pure noise, x (t) = n A {t), with a zero mean Gaussian probability 
distribution. The variance of this distribution is 

n M U)n B {f>) = l -P {A \f)5{f - f)5 AB . (16) 

By substituting for x A from Eq. (J2J) (with Z7 l GW and V^ GW set to zero there) in the above 
equation and using the covariances of the noise components defined in Eqs (jnj) and (|T5j) 
one finds [T4*] : 

3 

p {A \f) = E [d 2 ^ + r ?\ 2 + i 2 ^ A + ^| 2 )P proof 

i=l 

+ (\2p A \ 2 + \2q A \ 2 )P^] (17) 

The resulting noise spectra for each pseudo-detector is jS]: 

p(D (/) = p(2) (/) = 8sin 2 (7r/jL){ [2 + cos(27r/L)] P shot 

+ [6 + 4 cos(2tt/L) + 2 cos(47r/L)]P proof } , 
p( 3 )(/) = [2 + 4cos(2vr/L)] [P shot + 4 sin 2 (7r/L)P proof ] . (18) 

Therefore, the data analysis challenge is to detect signals in this remaining noise. 



3. The Signal 

Since LISA will be orbiting the solar-system barycenter (SSB), it is convenient to 
introduce a reference frame centered at the SSB. As shown in Fig. we define the 
x and y axes of this SSB frame to lie on the ecliptic, and the z axis to be normal to 
it and pointing towards the north ecliptic pole. The x axis points towards the vernal 
equinox. We take the GW source to be located in the direction given by the vector 
w. A gravitational wave from this source will arrive at the SSB origin traveling along 
— w. The sky position {6, <ft} defines the Cartesian components of the propagation 
direction, i.e., w = (sin 9 cos 0, sin 6 sin 0, cos 6). Thus, the sky-position angles are 



Optimal statistic for detecting inspirals with LISA. 



9 



equivalently characterized by two of the three components of w, say, W\ = sin 9 cos i 
and W2 = sin 9 sin </>. Also, 9 k and <j) k define the plane transverse to w: 







dw 



4» = - 



1 dw 



(19) 



89 sin 9 d(f) 

The perturbation created by the wave at a spacetime location (t, r) is given by 

h kl (t,r) = h + (t-w-r)(9 k 9 l -(f) k (j) l ) + h >< (t-w-r)(9 k (f) k + 9i(j)i) , (20) 

where h + (t) and h x (t) are the two GW polarizations and r is the position vector of the 
spatial location of the perturbation in the SSB frame. In the time domain, the strain 
induced along the ith arm is 

hi{t) = /i«(t)njn{ = h i+ (t)£ i+ (w 1 ,w 2 ) + h ix (t)^ ix (w 1 ,w 2 ) , (21) 

where we used the Einstein summation convention over the repeated indices k and I. 
Above, 

& + = (0 • n,) 2 -(</>• n^) 2 , &x =2(0-n i )(0-n i ) (22) 
are the beam-pattern functions of the ith arm for the two polarizations. 













* / A \ 

n/ \ \ 




' CCHpllG 


\ \ 


SSB 






• r ~~ip line of noiaV^^^^ 



Figure 2. The solar-system barycentric (SSB) frame, denoted by the (x, y, z) axes. 
The angles (9, 4>) specify the sky position, w of a GW source. The axes of the wave- 
frame associated with a gravitational wave emanating from that source are labeled by 
the unit-normal vectors (ic w ,y w ,z w ). The orientation of the wave frame relative to 
the SSB frame is given by the Euler angles, (</> — 7r/2, it — 9, ifi), where ip is the angle 
between the line of nodes and explicitized further on the left panel in Fig. |3J. 

The celestial longitude is drawn as a dashed arc passing through the origin of the wave 
frame. 

An impinging GW causes a change in the light-travel time along an arm that 
can be calculated by solving the null geodesic equation in the corresponding perturbed 
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spacetime. This in turn causes a time- varying Doppler shift, which clearly depends on 
the difference between the GW strains at the two space-crafts at the end of the arm. 
One also expects this shift to be dependent on the position of the source relative to the 
arm, w • fij. Thus, the GW contribution to the data stream Vi is given by 

-1 



Ht 



w • r m ; 



h(t 



w 



2(1 - w ■ n, 

where is the position vector of the ith craft in the LISA frame. 



Fi-l - L)\ 



(23) 




Elnecf notJe-4 




Figure 3. The left panel shows the wave frame and the right panel shows the 
orientation of the compact binary's orbit. The orientation of the wave frame is such 
that its z w axis points towards the origin of the SSB frame (as shown in Fig. |2J) and its 
x m axis lies along the semi-major axis of the wave's polarization ellipse. The tangent, 
1, to the longitude at the source's sky-position is perpendicular to the line of nodes, 
which lies on the ecliptic. The two form a plane on which x„, lies making an angle ip 
with the line of nodes. Therefore, (3 = n/2 — ip. In the right panel, COM is the center 
of mass of the binary. The normal, 6, to its orbit is along the binary's orbital angular 
momentum vector and has an inclination of angle e (introduced in the main text) with 
the line of sight, w. The projection of 6 on the celestial sphere is 6 p , which makes an 
angle (3 with respect to 1. The angle between 6 and 6 p is t. These two angles,/? and t, 
completely specify the orientation of the orbit's normal. Thus, i = n/2 — e. 

We now consider the effect of a signal from a non-spinning compact binary, with 
member masses m\ and m,2, on V^ GW (t). The two polarization amplitudes for the ith 
arm are 

1 + cos 2 e 



h i+ (t) = H(n i ) 
h ix (t) = H(Qi) 



2 

1 + cos 2 e 



cos 2ijj cos($i(t) + 5 C ) + cos e sin 2^ sin($ i (t) + S c ) 



(24,) 



sm2ip cos($j(t) + d~ c ) + cose cos 2-0 sin($i(t) +S C )\'25) 



2 

where <&i(t) is the phase received at that arm, {ip, e} are the polarization and inclination 
angles of the binary source (as explained in Fig. ^ and S c is the initial phase at the 
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origin of the LISA frame. For binaries that include a white-dwarf, the signal will be 
essentially monochromatic, with = flit. Above, H(Qi) is the signal amplitude at 
LISA defined as: 

2/3 



1.188 x 10" 



22 


" M 


5/3 


R 


-i 






1000 M _ 




lGpc_ 




27r x lmHz 



(26) 



where Ai = (m 1 m 2 ) 3 ^ 5 /(m 1 + m 2 ) 1//5 is the chirp mass, R is the distance to the source 
and 

a = ^o7* (27) 

is the Doppler shifted source frequency, Qq, at the ith arm owing to LISA's motion 
with respect to the solar-system barycenter. Note that 7$ is the sky-position dependent 
Doppler factor, 



li 



(28) 



where r\ i is the velocity of the geometrical center of the ith arm. Substituting the GW 
strain deduced from Eqs. (j^lj) and (}2*T]) in the V GW (t) expression (|2*3*j) . mainfests its 
dependence on the source parameters. One can similarly obtain Up w (t) as a function of 
these parameters. We first discuss the case of the monochromatic signal. The extension 
to slightly chirping sources is straightforward and will be presented subsequently. 

If the laser-frequency and optical-bench motion noises were at the level of the other 
noises, one could hunt for GW signals in the Ui and Vt data streams. But as we saw in 
the last section, this is not the case and, hence, one has to search for them in the pseudo- 
detectors x A {t) in which the contributions of these noises stand canceled. Working with 
the x A (t), however, makes the data analysis formulation a little non-trivial since one 
has to contend with time-delaying appropriately the six inter-craft data streams that 
can potentially harbor a GW signal, h A {t). Implementing this is somewhat easier in the 
Fourier domain. Thus, the algorithm we follow in the rest of this section is as follows. 
We begin by first computing the Fourier transform of the data streams V t GW (and U GW ). 
These will then be time-delayed and combined to calculate the frequency components 
h (u) (where ui = 2nf) of the GW signal in the Ath pseudo-detector, along the lines of 
Eq. ()13j) . Its inverse Fourier transform will finally yield h A (t), which is the quantity we 
aim to search for in the x A (t). In the process, we get an explicit expression for h A (ui), 
which is useful since the implementation of a search is faster in the Fourier domain, 
where one can avail of the existing Fast Fourier transform algorithms. 

We begin by defining two new functions of {ip, e} that appear naturally in the 
Fourier transforms, V^ GW (cu) and U\ ; (a;), of the six streams: 

l x = -i(T 2 2 (^e,0)-T 2 - 2 (^,e,0)) , 

/+= (T 2 2 (^,e,0) + 77 2 (^,e,0)) , (29) 
where T 2 ±2 are Gel'fand functions 



rp ±2 
1 2 



e, 0) = - (1 ± cos e) 2 exp (=Fi2^) • 



(30) 
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To find V? w (u), we use the GW strain Eqs. and (HQ in the expression for V? w {t) 
in Eq. (}2*3*|) . Taking the Fourier transform of the result gives 



V™(lu) = H(n )T 



sinc((w - Vti)T)e~' l5c x 

(/;(</>, 6)F* +i;+ (^) + T x e)^ +i;x (^)) 
+ sinc((o; + n i )T)e i5c x 

(Z + (V,e)F ym;+ (a)+ix(V' ) e)F yi+1 . x (fi i )) . (31) 

The orientation of the ith arm resides in the functions, 

Fv i+r ,+,x = -i&i6;+, x sinc (^L^r) e iT " l ^ L! , 

F^.+.x = ifoi^xsinc , (32) 

where 

^^4^, Ti^ifi-^y ^= (i± f ^ , (33) 



2 V 

are all real quantities. The Fu { ;+,x similarly determine the fractional frequency shift 
[/ Gw . Note that the depend ence on the angles {ip, e} has been separated out in the 
form of /+.x- We will exploit this separation of variables in the next section to eliminate 
the computational cost in searching over the {ip, e} angles for a GW signal. 

The Fourier components of the GW strain h A {u) can now be found by combining 
the above V^ w (u) (and U^ w (cu)) via the Fourier analogue of Eq. (JT3J). Thus, 

3 r 

e -i(5 c+CT f) sinc _ T ) T p* D M e -i n ciiLi 

, (34) 



h A (u) = iH(n )Tj2 



+ e i(5c+ff * Wc ((u + T) T 2 " D A e 1T ^ iLi 

where there is an implicit sum over p = ±2. Also, we define 

DLi = bj\Mf\ (C. I iC, ) , (35) 

where 

M A = q A lS mc (QiLikt) - p? +1 smc (fi^^r) (36) 
and of = arg(M^). 

The time-domain expression of the GW strain in the pseudo-detector A is obtained 
by taking the Fourier transform of h (uj), and is found to be 

3 

h A (t)=H(n )J2^[e- i5 Ef* S,(t)} , (37) 

where 5 = 5 C + ir/2, 

S A (t) = e in ^- L ^ )+ia t/ g A and E A = g A T 2 p D A . (38) 
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Above, is a normalization constant such that 



A A ^ D J W > 

(Sf,Sf) (A) =4 j[ rf/^y^ = l , (39) 



which implies that for an observation duration (i.e., signal integration time) T, 



4 



IT 



lV2 



(40) 



When considering a slightly chirping source (i.e., when Cl Q <C VLq/T) one can expand 
the signal phase as 

§ l {t) = a i t + Ui i t 2 . (4i) 

In that case, the only modification to h A {t) occurs in the Sf term: 

Sf(t) = ^i( t - L 3 T i)+ i (k* l i#+°£) /gA ( (42) 

which defines the time-domain template for chirping compact-object binaries. Note that 
apart from the normalization constant, the template SHt) is a pure phase term. We 
will find this useful when deducing the matched filter in the following section. 

Equations (j2ZJ) and (J2%|) show that flj is related to the sky position and the intrinsic 
chirp rate Cl . In order to ensure that the phase evolution of a template models that of 
a signal well, it is important to allow for non-zero m searches of binaries involving 
masses higher than those of white dwarfs, such as in searches of neutron-star binaries 
and binaries involving a neutron star and a black hole. This is because the gravitational 
radiation reaction on these sources is stronger than those involving white-dwarfs. The 
post 1 -Newtonian waveform reveals that for small chirp masses and source frequencies, 
the waveform phase can be expanded as in Eq. (jUJ) , with 

Thus, determining Qq is significant since, together with Qq, it determines the binary 
chirp mass A4. And as shown by Eq. ()26|) . additional knowledge of the amplitude will 
then help in estimating the distance to the binary. ^5] 

To summarize, the GW strain in pseudo-detector A is given by Eq. (|3*4*|) and is 
determined by eight independent parameters, {R,5,Qo,Clo,ip,e,9,4)}. To search for a 
signal we must devise a strategy to seek these strains for a range of parameter values 
accessible to LISA's pseudo-detectors. This is what we deal with in the next section. 



4. The Optimal Statistic 

Given three independent pseudo-detectors, x , we now ask what is the optimal detection 
statistic to look for GW signals, h A , in them. In the absence of any prior probabilities 
and costs, the optimal detection strategy is the one that minimizes the rate of false 
dismissals for a given rate of false alarms. This is termed as the Neyman-Pearson 
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criterion. Under this criterion, and for zero-mean Gaussian noise, the detection statistic 
is the likelihood ratio, A, defined as 

A ,x A ) {A) -l{h A ,h A ) (A) ) . (II) 



A=l ^ 



where the first term is the sum of the cross-correlations of the expected signal, h A , with 
the respective data, x A , over all pseudo-detectors. The cross-correlation for pseudo- 
detector A is given by 

{h A , x A ) {A) ee 4K jT df hA J^ f) , (45) 

where 3ft (X) denotes the real part of a complex number X. The second term in Eq. (}4*4*j) 
is an overall normalization that is independent of the data. Substituting for h A from 
Eq. (J3ZD, we get: 

A^E^t^^O*] > ( 46 ) 

A =l A =l i=l 

where 

C A =(S A ,x A ) (A) . (47) 

The double summation in Eq. (J46|) can be replaced with the single sum over a new 
index k, 

3 3 9 
A=l i=l k=l 

where A = Ceiling(A;/3) = [§] and i = Mod(fc-l, 3) + l = (fc-l)%3 + l.|2H We simplify 

the above expressions further by consistently identifying Yn^x^xs+i with Y k , which form 
the components of a 9-dimensional vector Y. Thus, Y k=1 = Y^ 1 , Y k=2 = Y^ 1 , 
yfc=3 = Y^f 1 , Y k=4 = Y { A = 2 , and so on. We use this algorithm to map the D A , E A , 
and C A to the components, D k , E k , and C k of 9-dimensional vectors D, E, and C, 
respectively. 

In the new notation, the cross-correlation statistic becomes 

3 9 

J2(h A , x A ) (A) = & b~ iS E* k C k ] = 3ft [e- iS E • C] (49) 



A=l k=l 



and 



J2(h A , h A ) {A) = H 2 (n ) E* k E k = H 2 (n ) || E || 2 ee k 2 , (50) 

A=l k=l 

where || Y || denotes the norm of vector Y. Therefore, k 2 is a measure of the signal 
power accessible to LISA. It is usually less than the peak power H 2 (fl ) owing to LISA's 
non-optimal orientation, E, to a given source. The relative sensitivities of the three 
pseudo-detectors to different sky positions can now be studied by plotting the analogue 
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Figure 4. Sensitivities za = E A * Ef / (gf) 2 of the three pseudo-detectors as 

functions of sky positions, {8, </>}, in radians. Note that the Doppler shift in the source 
frequency has been ignored here; so gf is the same for all i, i.e., gf = g . These 
plots have been evaluated for CIq = 1 mHz at the "initial" orbital position of LISA 
labeled t — 0. As illustrated above, direction of maximum sensitivity varies from one 
pseudo-detector to another. 



of k 2 for each individual pseudo-detector, as shown in Fig. HJ This figure, plotted for 
Qq = 1 mHz, verifies the fact found in Ref. |B] that the third pseudo-detector has a 
much smaller sensitivity than the first two. It also, shows that the peak sensitivity of 
the pseudo-detector labeled as A = 1 is the best of the three. However, at any given 
location on LISA's orbit, there are sky positions to which pseudo-detector 2 has the best 
sensitivity. We also plot in Figs. El and El the sensitivities of pseudo-detector 1 and the 
optimal combination of all the pseudo-detectors, for three different locations on LISA's 
orbit. These figures show that for all sky positions the optimal-sensitivity is better than 
the usually best pseudo-detector, labeled A = 1. 

The likelihood ratio now takes the following form: 

lnA = K X>V A >(A)-^ 2 , (51) 



A= 



where 



h A Lj= - | Q ■ C | 2 , (53) 



is the normalized counterpart of h A , such that Xm=i(^ A > ^ A )(A) — 1- The likelihood 
ratio can be maximized over k and 5 to yield 

1 
2 

where the hat on a parameter denotes its value at which In A stands maximized with 
respect to that parameter and Q is the normalized orientation vector, Q = E/ || E ||, 
such that || Q ||= 1. Also, we find k = J2A=i(h A ,x A )(A) and 5 = arg(C • Q). 

To maximize the statistic in Eq. (p)3*|) with respect to {ip,e}, note that Q can be 
expressed in terms of its components as follows: 

Q = Q +2 D +2 + Q- 2 D_ 2 , (54) 
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where 



±2 ,Jj^ and (55) 



Above, we have used the fact that || D + || = || D || = || D ||. The statistic in Eq. (JHHj) 
depends on {ip, e} solely through Q. Therefore, it stands maximized with respect to 
those parameters when Q gets aligned with C. The fact that this alignment is physically 
realizable was shown in Ref. ESI- The maximized statistic is 



= i || C„ || a > (56) 

where Cn is the projection of C on a 2-dimensional complex space, TC, spanned by 
{D +2 ,D_ 2 }. Since it is always possible to choose a pair of real basis vectors to define 
this two-dimensional space, we take these vectors to be 

6 ± = (di ± d 2 ) I || di ± d 2 || , (57) 



where di = 9?(D +2 ) and d 2 = S(D +2 ). Thus, we may define the network search statistic 
as 

A =|| C n \\ 2 = \C+\ 2 + |C-| 2 = (c+) 2 + (c+ 2 ) 2 + (c - ) 2 + (c; /2 ) 2 , (58) 
where = o 1 * 1 • C = + ic^, 2 . The maximizing values of the two parameters are 

V> = arg(ro)/4 and e = cos -1 (1 - A /|ro|)/(l + a/M) , where w = C +2 /C~ 2 . 

The above statistic is still a function of the sky position through the parameters 
{wi,w 2 }. Since it is not possible to maximize the statistic over these parameters 
analytically, one must resort to a numerical maximization scheme as described in the 
following section. By comparing the values of the statistic for each pixel in the sky 
with a threshold value, A , a decision on the presence or absence of a signal in the data 
can be made. The threshold itself is determined by the false-alarm probability that one 
can afford. Note that in the absence of a signal, A is a random variable that has a \ 2 
probability distribution, 

j9 (A) = ^exp(-A/2) , (59) 

with 4 degrees of freedom [18 . This is because each of the Cq and cr, 2 is a Gaussian 
random variable with a zero mean and a unit variance. The false-alarm probability is 

Qo = ^°°R,(A)d(A) = (l + exp (-Ao/2) . (60) 

In the presence of a signal, the probability distribution of A is non-central x 2 ; 

Pl (A) = ^exp(-^)/ 1 ( K VA) , (61) 

with the non-centrality parameter as k 2 , which is a measure of the signal power j!8j . 
Above, I\ is the modified Bessel function. 
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Figure 5. Sensitivity zi, as defined in Fig. 0J evaluated at fio = 1 mHz for three 
different angular orbital positions (in radians), {0, ^, with respect to the "initial" 
location denoted in Fig. 0] The left plot above is identical to the left plot in Fig. 
0] since it corresponds to the same pseudo-detector and orbital location. Note that 
the sky positions corresponding to the sensitivity maxima vary from one location to 
another on LISA's orbit. 




Figure 6. Network sensitivity z op t = J2t=i E%.E k / (g 1 ) 2 evaluated at fio = 1 mHz for 
the same orbital positions that appear in Fig. |5] It is manifest that pseudo-detector 3 
makes negligible contribution to the z op t at this frequency. Note that g 1 — g 2 . At any 
given sky position, the optimal statistic has better sensitivity than any za- 



5. Template Spacing and Computational Costs 

The A statistic in Eq. must be maximized over the remaining intrinsic parameters, 
namely, •& = {Qq, fio, w\, ^2}- As a first calculation, we will set fio = and will focus on 
the number of templates required to scan the space of {Oo, w±, W2} for this case. As noted 
above, the ensuing template spacings will still be relevant to a large number of compact 
binaries that involve white-dwarfs. The corrections arising to these spacings for non- 
zero Cl will be studied elsewhere. The maximization can then be achieved numerically 
using a discrete template bank over this three-dimensional parameter space. The drop 
in the value of the statistic and, therefore, the signal-to-noise ratio (SNR) that one can 
afford determines how coarsely one can space the templates. In practice, there are limits 
posed by the available computational resources on how fine the spacing can be. The 
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loss in SNR is related to the template spacing through the ambiguity function. [20] 

The ambiguity function corresponding to the A statistic is derived from it by 
replacing the data x A there by a signal We distinguish the signal parameter 

values from those of the template by denoting the former with a prime. The parameter 
values of a template used in a search may not be the same as those of a signal hiding in 
the data. Let the parameter mismatch be A# = (-&' — Then the ambiguity function 
is a real quantity expressed as 

m(#,Atf)= Pl !Q' k Q' l *e mh (62) 

where Q' k depends only on the signal parameters and 

= e i * fe e Wi sinc(fifcT - fi fc T)sinc(^T - ftyT), (63) 

and p k l = o +l o k + o~ l o k is an amplitude factor. It is important to note that as Ai9 — > 0, 
one has m(i9, Ai?) — ► 1, which is the maximum value it can attain. The correlation 
phase, ^k, is defined as 

= L{VL k T k - Q k r k ) = Qk-^^Awx + Vl k ^=Aw 2 - r k AQ k . (64) 

The drop in the value of m(i9, Ai9) caused by non-zero, but small A-$ M , can be 
ascertained by Taylor expanding it about the maximum at A$ = 0. [221 123] The first 
order term is zero since by definition the statistic has a maximum when the template 
parameters match the signal parameters. Thus, 

1 - m(0, Atf) ~ ^j a pAd a Ad 13 , (65) 

7 Q/ 3 is determined from the second order term in that expansion: 

1 /<9 2 m(tf,Atf) 



1af3 



2 VdAtf Q <9A^ 



(66) 

Al9=0 



It defines the metric on the 3-dimensional parameter space. 

The computational cost for the search can be reduced by taking advantage of the 
Fast Fourier Transform algorithms |17. and computing the cross-correlation components, 
C, in the Fourier domain. This defines the strategy for searching for the source 
frequency, To search for the remaining parameters, {wi,W2}, one must design a 
bank of "templates" with values of these sky positions spaced such that the loss in 
SNR is never more than the desired fraction, say, [l. To find the metric, on the 
two-dimensional space V, spanned by {wi,^}, we project 7 Qj a orthogonal to Qq, 

9ij = lij , (67) 

7oo 

where i and j span only the {wi,Wz} space, and the index denotes the flo ax is- The 
volume of a V is then given by 

V= f y/det\\g ab \\d p ti , (68) 
Jv 
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where P = 2 is the dimensionality of the space. The number density of templates, 
pp(p), is determined as a function of p to be [25] : 

M = (l]jf) ■ (69) 

Setting the fractional SNR loss p = 3% yields a template density of p 2 (0.03) = 16.6. 
Therefore, the total number of templates is just the overall parameter volume times the 
template density, i.e., iV temp i ates = V x p P (p). 

The parameter volume, obtained via the metric computation, turns out to be about 
5 for Q = 1 mHz. Considering that the template density per unit volume is only 16.6 
implies that the number of sky-position templates required for a search with 3% loss of 
SNR is about 80. The smallness of this number is not surprising, given how slowly the 
ambiguity function varies as a function of A6 and A<j), as shown in Fig. [7| 

As is manifest from Eq. (|62jl . in principle, this variation can arise from either the 
time delays in Q^i or the weights p k l . However, for wavelengths much larger than the 
LISA arm length, the ability to discern between different sky positions through the time 
delays in B^ is negligible. The main contribution to m(i?, Ai9), therefore, arises from 
the P)}. For detailed studies of the angular resolution achievable by LISA, we refer the 
reader to Refs. jMl23]. 




Figure 7. The ambiguity function m plotted as a function of the template parameters, 
{9,4>}, for three different source sky positions (in radians), namely, {9',<j)'} = 
{7r/2,47r/3}, {7r/6,7r/3}, and {n/2,Tr}. 

To get a handle on the computational costs associated with a search of this nature, 
one must determine the overall number of sampling points in a single data train. For 
our analysis, we have chosen a sampling rate of 1 Hz. Therefore, the number of sampling 
points is just the length of the data train, T. The number of floating point operations 
associated with an FFT with T sampling points is: 

iV fl _ opts = Tlog 2 T. (70) 

Considering that this must be done for each template, the resulting expression for 
the total number of floating point operations for an arbitrary template bank is just 
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-^templates x ^V"ft-opts- However, the real quantity of interest is the the number of operations 
per second, 

C flops = iVtcmplates T X Nh ~ opts = AW plates log 2 T • (71) 

Therefore, the computational costs of implementing a real-time search with about 
80 templates on a year's worth of data is trivial since the number of flops scales 
logarithmically with the integration time. 



6. Discussion 



In this paper, we developed an optimal method for detecting slightly chirping compact- 
binary inspiral signals in the LISA data. We also studied the geometrical properties 
and sensitivities of the three noise-independent pseudo-detectors or data combinations 
of LISA. Following the earlier work on TDI data combinations, it was found fTJ EI] 
that the rotational motion of LISA and the time-variation of LISA's arm-lengths 
would prevent the noise contribution of the laser-frequency fluctuations from being 
mitigated to the level of the secondary noises. In order to tackle this problem, 
second generation pseudo-detectors were introduced as simple differences of their first 
generation counterparts, appropriately time-shifted: 

x A (t) = x A (t) - x A (t — Li — L 2 — L 3 ) ~ x A (t) - x A (t - 3L) . (72) 

The analysis presented above can be easily extended for detecting inspiral signals in the 
second generation pseudo-detectors, x A (t), by implementing the next two steps: First, 
by following the derivation in Sec. |2] it can be verified that the noise PSDs of the new 
detectors are given in terms of the old ones by: 

p (i,2,3)( / ) = 4sin 2 (3vr/L)p( 1 ' 2 ' 3 '(/) . (73) 

Second, since a GW contribution to x A {t) will also get differenced as in Eq. (JZ2J), the 
template that must be matched against x A {t) should itself be modified accordingly: 

Sf(t) = Sf(t) - Sf(t -L,-L 2 - L 3 ) ~ Sf(t) - S A (t - 3L) . (74) 

None of the formal analysis presented in this paper is affected by this change. For 
example, the analytic maximization of the likelihood ratio over four source parameters 
(namely, the signal strength, the initial phase, the polarization angle and the orbital- 
inclination angle), and the concomitant computational gain achieved in the process, still 
hold so long as the matched-filter output in Eq. (jl7|) is redefined to be 

C A =(Sf,x A ) (A) . (75) 

Also, the sensitivity plots in Figs. El and El for the first generation detector A are 
the same as those of its second generation counterpart, i.e., detector A. This is because 
the (geometric) sensitivity, za, depends on the orientation of a detector relative to the 
source and is independent of the noise PSDs. Similarly, the formal expressions for the 
ambiguity function in the last section above remain unchanged. But, since this function 
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depends on the noise PSDs, its numerical value is affected by the change in Eq. (J73J). 
Nevertheless, we have found that this change has negligible effect on the number of 
templates deduced in Sec. El 

A useful by-product of our analysis is that it yields the maximum-likelihood 
estimates of the initial phase, the polarization angle, and the angle of inclination (in 
addition to, of course, the signal amplitude). To complete the parameter estimation 
problem, however, one needs to derive the errors associated with them as well as in the 
parameters that will be searched for numerically, viz., {Oq, Qq, 6, 0}. That problem will 
be addressed elsewhere. 
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